A Micro-Thermodynamic Formalism 



Hans Henrik Rugh 
Department of Mathematics, University of Cergy-Pontoise, 

France 95302 

February 8, 2008 



Abstract 

We consider the micro-canonical ensemble of a classical Hamiltonian dynamical system, the Hamil- 
tonian being parameter dependent and in the possible presence of other first integrals. We describe a 
thermodynamic formalism in which a 1st law of thermodynamics, or fundamental relation, is based 
upon the bulk-entropy, Sn- Under an ergodic hypothesis, Sn is shown to be an adiabatic invari- 
ant. Expressions for derivatives and thermodynamic relations are derived within the micro-canonical 
ensemble itself. 

Equilibrium properties of an isolated Hamiltonian dynamical system with many, say lO^**, degrees of 
freedom is probably best described using a thermodynamic formalism for a canonical ensemble at a fixed 
temperature, even though this means introducing fluctuations in an otherwise conserved quantity, the 
energy. For a more moderate number of degrees of freedom, say 10'^ — 10®, numerical simulations become 
feasible, and it is desirable to obtain a description in terms of the micro-canonical ensemble itself where 
the values of the first integrals are fixed quantities. In such an approach, geometrical properties of the 
level surfaces reflect thermodynamic relations and, by invoking the ergodic hypothesis, also dynamical 
properties of the underlying system. In particular, when energy is the only first integral, measurements 
may then be done by time averaging (cf. ||l^, |ll[). The purpose of this article is to develop such a micro- 
thermodynamic formalism further, taking into account parameter dependency and the presence of other 
first integrals. Within this framework we will also (section ^) discuss a natural formulation of a 1st law of 
thermodynamics, or fundamental relation, based upon the bulk-entropy, Sq. We refer to Abraham and 
Marsden Q as well as Landau and Lifshitz ||^ for a general introduction to thermodynamic ensembles, 
to Lebowitz et al. for an illustrative example of some differences between the ensembles and to Evans 
and Morriss H for practical calculations carried out in the micro canonical ensemble. We also refer to 
Jepps et al. 0^ where the presence of other (approximative) first integrals is of relevance and to Otter 1^ 
who studied reaction events using mixed ensemble averages. 

1 Micro-canonical ensembles 



For simplicity, we consider a Euclidean phase space, fl — R , d > 1, and a Hamiltonian function, 
H : Q ^ R, bounded from below and of sufficient rapid growth at infinity. The dynamics preserves the 
Liouville measure, here the Lebesgue measure, m — d'^cflp. There may be other first integrals, denoted 
F — Fi,. . . ,Fm, m > Note that in section |^ we shall write Fq — H for the Hamiltonian which is 
then considered at the same footing as the other first integrals. All first integrals are assumed to be in 
involution. We also assume that all functions are known analytically and that the Hamiltonian depends 
smoothly on some external real parameters, denoted A = Ai, . . . , A„. By contrast we do not allow the 
other first integrals to depend on A. This is for technical reasons (cf. below), though in some cases such 
a condition could be relaxed. For fixed values of parameters, A = A, of first integrals, F — I, and of the 
energy, Hx — E, the subspace, A = A[E, I, A] = £ i?^'^ : Hx{S,) — E, F(^) — I}, is invariant under the 
dynamics of Hx- We will assume that values are chosen so that the differentials, dH, dFi,. . . ,dFm, are all 
independent on A[E, I, A]. This in particular implies that A is a smooth co-dimension m+l sub-manifold 
of our Euclidean space. 

In the literature one will find (at least) two definitions (denoted bulk and surface) of a micro-canonical 



entropy and temperature. It turns out that we shall need both. Thus we define 



where Q denotes the Heaviside function and 



j mQ{E-H^)5{I-F) (1) 
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The bulk- and the surface-temperature are then given by: 

Tn dE' dE' ^ ' 

Derivatives with respect to other first integrals are considered in Section (||). We also have generalized 
bulk- and surface-pressures: 

i rr 9Sn j dSf^ . , , 

Pn^Tn-^, Pf^^^f^-Qx^' « = 71. (4) 

Taking an average in the micro-canonical ensemble will here mean taking the surface-average, i.e. 

J md [E-Hx,I-F) 

In practice, {4>)^, is often calculated by time-averaging (assuming ergodicity), thus giving a dynamical 
preference to the surface-average relative to other ensemble-averages. 

The calculation of either of the two entropies may be difficult or even impossible when the number of 
degrees of freedom in the system is large. On the other hand, the associated temperatures and generalized 
pressures may be probed using time- averaging (cf. Section A bulk-pressure, Pq, may be calculated as 
follows: 



Pii = -{^),^~{^\EJA)- (6) 



To see this we note that the derivative of a Heavyside-function yields a delta-function. It follows that 
l/Tn = dSn/dE = e^'^/e'^'', and therefore, 

_^ dSn JmSjE-HxJ-F)^ 
^""^''9A7" JmS{E~Hx,I~F) ' 
By the very definition, there is always a 1st law of thermodynamics for the bulk-entropy: 

TndSn = dE + Y,PlidK. (8) 
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The same kind of relation is, of course, valid for the surface-entropy but, as we shall see, adiabatic 
invariance clearly gives a preference to the version involving the bulk-entropy: 



2 Adiabatic invariance 

An adiabatic process is characterized by a slow and smooth change in parameters during which 'thermal 
equilibrium' is maintained. For example, slowly moving a piston of a cylinder containing a gas of particles. 
In order to to avoid what is known as parametric resonance, we fix a smooth path in parameter space, 
s G [0, 1] 1-^ c(s) € i?", and traverse this path in rescaled time, i.e. t G [0, A] i-^ X{t) = c{t/A), for some 
A > 0. A physical trajectory, ^{t), is evolved from an (almost) arbitrary point, ^(0) £ A[E, I, A], using the 
time- and A— dependent Hamiltonian, ^fA(t) ■ The goal is to determine the micro-canonical state, i.e. the 
values of {E, I), of the system at time t — A. The process is called adiabatic if this state has a well-defined 
limit as A — > oo, the limit depending on the initial state {E,I,X) and the path, but not on the choice 
of initial point, cf. Arnold ||^, section 52, or ||l^. Here again, we do not permit the other first integrals 
to depend on A. Being time-independent and commuting with the Hamiltonian, they therefore remain 



constant under the time-involution. The energy, E{t) ~ Hx(t){S.{t)), being time-dependent through A, is 
in general not constant in time. Instead we get by Hamilton's equations: dE/dt = dH\Jd\i ■ dXi/dt, 
and the 'adiabatic' problem is then to see if E{A) has a well-defined limit as A ^ oo. 

In order to attain the adiabatic limit, A should be very large, in particular larger than the time-scale, 
Torg, over which ergodic averaging takes place. We may then look at the time-average of an observable 
over an intermediate time-scale, Tcrg << St << A. This time-scale is short enough so that A does not 
change significantly but long enough so that a dynamical average may be replaced by a micro canonical 
average. Under this 'Adiabatic Ergodic Hypothesis' we get for the energy-change: 

Over the time-scale, 6t, the A-derivative is almost constant (because of the scaling with A) and may 
therefore be taken outside the average: 

i 

an expression which is correct to order Tcrg/A. Taking the A ^ cxd Hmit we get an identity between 
differentials, 

dE = Y,{^\EJ,\) dK (11) 

i 

valid precisely in the adiabatic limit. By (|^) and (H), we see that the bulk-entropy, Sn, is indeed an 
adiabatic invariant. In l-fl dimensions the bulk-entropy is just the action integral and the phenomena is 
well-known (adiabatic invariance of the action, cf. Arnold [^). In higher dimensions a similar result was 
obtained by Kazuga though in a different context. 

Now, the bulk-entropy is a fortiori strictly increasing as a function of the energy (at given parameter- 
and first integral- values) . One particular gratifying, though non-trivial (!), consequence is that by travers- 
ing a loop adiabatically the energy must return to its original value. 

It is easy to put the Adiabatic Ergodic Hypothesis into a more rigorous form (in terms of decay 
of correlation functions). It is, however, virtually impossible to check analytically if such a condition 
really holds in a given situation. One serious problem is that critical slowing down (meaning that Torg 
diverges) occurs if one encounters an additional first integral along the traversed path. On the other hand, 
singularities in the energy surface are likely to pose problems only in low dimensions (notably 1+1). 

Should one allow other first integrals to depend on the parameters, the bulk-entropy, as we have 
defined it, is in general no longer an adiabatic invariant. We have not been successful in finding a good 
replacement for the bulk-entropy (the reader is encouraged to try for himself) and consequently not 
allowed such a parameter dependence. 



3 Surface - derivatives 

Close to equilibrium we may express response functions as derivatives of averages with respect to parame- 
ters and values of the first integrals (in this section this includes the Hamiltonian). When the Hamiltonian 
is the only first integral and there is no parameter dependence, it was shown in |o[ |ll| (cf. also |, |]) that 
energy-derivatives of a micro-canonical average may themselves be calculated as averages. The general 
case turns out to be quite similar. 

Our assumption on the first integrals being independent means that in a neighborhood of the integral 
surface, A[E^ /, A], we may construct vector fields, Xq,. . . ,Xj„, for which 

dF^{X,) = 6^,, i,j ^0,...,m. (12) 

Geometrically, the vector field, Xi, is transversal to the Fi = li surface but parallel to the other surfaces, 
Fj = Ij, j ^ i. We make the observation that if f{x) is a (suitably smooth) function and V is a vector field 
for which df{V) = 0, then the Lie-derivative, LvS{f{x)) — {V ■ V)6{f{x)) = 0, vanishes identically. This 
is clearly true if the delta-function were a smooth function and the claim then follows by approximation 
(a rigorous proof is quite lengthy). This applies to our Hamiltonian set-up, since for all i, j ~ 0, . . . ,m 



Hence, if (j) is any smooth function, 

y ™ • + • V))<5(/-F) = 0. (14) 

Taking the /^-derivative outside the integral and carrying the Lie-derivative out by partial integration we 
obtain 
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J m 6(1 - F) <j) = j m5{I-F) V-i<j)X,). (15) 



The identity, (gj^ -I- ■ ■^)S{I — F) = (recall that only the Hamiltonian depends on A), also 
implies: 

A_jTr,5{I-F)4> = -Jm 5(1 -F)\7. i^^Xo). (16) 

The surface entropy is the logarithm of a surface integral. Hence, when taking derivatives a normali- 
sation factor appears which precisely turns the derivative into a micro-canonical average. Thus, for the 
inverse 'generalized' surface-temperatures, = (with = l/Tf^), we get using ( p5[ ) and setting 

0=1, 

^ f m S(I -F) (j) 
^'^= )m8iI-F)^ HV.X.|i.,/,A), (17) 

Similarly for the 'generalized' pressures, using (16[), 



When taking a derivative of an average and writing 5 — 5{I — F)) we have 



qt'Wa'^ f — X r — X f — x~ (19) 

oli J mo jmd j mo 



which by the definition of averages and generalized temperatures reduces to 

_d_ 

We calculate in the same way. 



(0)^ = (V-(0X,))p-/3;,(<^)^s. (20) 



— {cf,), = -(V • (0^^o)), -pU^),. (21) 



These thermodynamic identities provide the natural generalizations of the results found in ||T^, ^ . In the 
above formulae the other first integrals may, in fact, be allowed to be parameter dependent (essentially 
because there is no Heavyside function in the above). The straight-forward derivation of formulae in this 
case is no more difficult and left to the reader. 

We also note that if the Liouville measure has the form dm = piO'^'^'^^ (fo^' instance in local coordinates 
on a symplectic manifold) , the only change in the above formulae is to replace the divergence of a vector 
field, V • by iV • (pF) (cf. Ilj). 

4 Bulk - derivatives 

Should one wish to take derivatives in the bulk-ensemble the procedure is slightly different. Restricting 
our attention to the energy derivative, suppose that F is a smooth vector field such that [compare with 

V-r = l and dF,{Y) = 0, I <i <m. (22) 

One may certainly find such Y when there are no other first integrals present. In the general case it is 
less clear because here we need the vector field to be defined 'smoothly' throughout the 'bulk'. Assuming, 
however, that we have found such a vector field we note that 



J m (V • r) e{E - H)6{I - F) (23) 



which by partial integration yields 

2^" = I m dHx{Y) 5{E - H,I - F) (24) 



e 

In particular, we obtain the following (well-known, when energy is the only first integral) formula for the 
bulk-temperature 

Tn = e^"-'^- = {dHxiY)\E,I,X). (25) 
For comparison we note that in the canonical ensemble where an integral of (f> looks as follows: 

m e-f^"^ S{I - F) 0, (26) 

the same calculation shows that 

Tcan = = {dHx{Y)\f3,I,X). (27) 

In the literature dH(Y) is therefore often used to 'define' the micro-canonical temperature (again mea- 
sured by time-averaging). With no other first integrals present one may take Y to be the canonical 
momenta divided by the number of degrees of freedom, i.e. Y = P/d. In that case, dH{Y) — i(iiJ(P) = 
l/d ^ (jiPi is the normalised reduced action (proportional to the kinetic energy, when quadratic in mo- 
menta). By time-averaging, we obtain a quantity proportional to the reduced action-integral. This is 
extremal under variations of the trajectory preserving the energy (Maupertius principle, Section 45]). 
One would thus expect a finite time-average of dH{P) to 'probe' a larger neighborhood of the trajec- 
tory than time-averages of other observables. This suggests a faster ergodic averaging, whence stronger 
numerical stability, when calculating time-averages of dH{P) relative to other observables of the same 
ensemble- variance. 

5 An example 

Consider an ensemble of N particles moving on a 3-dimensional torus, {R/Z)^, under the influence of 
pair-potentials. Putting things on a torus compactifies the configurational space but otherwise does not 
affect our results. The Hamiltonian, H = X]i3=o Pi Z^'™* + Si j U{xi — xj), is translational invariant, 
hence the total momenta, ptot — ^Pii provides 3 first integrals in addition to the Hamiltonian. If one 
here uses the normalised total kinetic energy as a measure of the temperature one runs into the following 
paradox: A configuration which minimizes the potential energy but has all particles moving at the same 
constant velocity, vq, is stationary. With no 'apparent activity' one should assign a temperature zero to 
this configuration but this is clearly not what the total kinetic energy does. 

One solution to this paradox suggests itself: Move to a frame where the average velocity vanishes. 
Such a procedure is, in fact, a natural consequence of our micro-canonical formalism. To get the pre- 
factors right we carry out the details: If P = (pi, . . . , p^v, 0, . . . , 0) denotes the full momentum vector in 
i?"^^ X {R/ZY^ and v = X^jPi/X^i'^i the average velocity then the reduced momentum vector and 
kinetic energy is given by Z = (pi - miv, . . . , pw - mArv, 0, . . . , 0) and K^ed = I]i(P* ~ TOiv)^/2mi, 
respectively. Straight-forward calculations give dHiZ) = 2K-ccd and dTptot{Z) = 0. The latter is precisely 
the requirement of 'parallelism' needed for both ( [l^ ) and ([2^). Furthermore, we have V-Z — 3(iV— 1), 
and \/ -{Z / dH{Z)) = (3(7V— 1) — 2)/(2ifiod)- After the appropriate normalization, formulae ( [l7| ) and ( |25| ) 
yield the surface- and bulk-temperatures: 

-1 = (^(^-^)A and T. = (4^),. (28) 

In both cases the conservation of 3 momenta results in a subtraction of 3 degrees of freedom as compared 
to the unconstrained case (pleasing on physical grounds). It is interesting also to note that the above 
formulae involving the reduced kinetic energy are valid irrespective of the total momentum being preserved 
or not ! When the total momentum is only approximately preserved one could use the above construction 
to define an approximate 'instantaneous' and possibly 'local' temperature. In [Q, such a problem was 
considered and the solution suggested was to choose a vector field, Z , depending on the configurational 
coordinates only. This automatically ensures 'parallelism', i.e. d'ptot{Z) = 0, and leads to what is often 



denoted the 'configurational temperature'. On the other hand, for numerical reasons the reduced kinetic 
energy may be preferable. 

The case of a conserved angular momentum (an experiment involving an axial symmetry) follows the 
same straight-forward procedure and is left to the reader. 
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